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We present a general technique to compute how the energy ol a configuration varies as a lunction 
of its overlap with the ground state in the case of optimization problems. Our approach is based 
on a generalization of the cavity method to a system interacting with its ground state. With this 
technique we study the random matching problem as well as the mean field diluted spin glass. As 
a byproduct of this approach we calculate the de Almeida- Thouless transition line of the spin glass 
on a fixed connectivity random graph. 
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I. INTRODUCTION 



The study of the ground state properties of disordered systems reveals deep connections with the field of random 
combinatorial optimization: in fact combinatorial optimization problems can be stated in terms of zero-temperature 
statistical properties of some disordered system Hamiltonian. A closely related issue is the computational complexity 
of a given problem that can be stated in terms of the typical amount of time {e.g. cpu-cycles) that is needed to 
solve the problem as a function of its size p). Easy/Hard thresholds have been observed in different combinatorial 
optimization problems @, 0, 0, 01 j ancl a currently highly debated issue is the interplay between the onset of the phase 
transition and the slowing down of local search algorithms A better understanding of the organization of the 
lowest energy configurations is then a promising research program not only to relate the computational complexity of 
a given problem to its statistical mechanics counterpart, but also as a guide for the implementation of more efficient 
algorithms. 

Recently the problem of diluted mean field spin glasses received a renewed interest by means of the cavity approach 
0, U- Interestingly this method turned out to be a powerful tool also to deal with combinatorial optimization 
problems. In this paper we study the problem of the organization of the lowest energy configurations within this 
framework. Our aim is to explore the way the energy of a configuration varies as a function of its overlap with the 
ground state. The naive procedure consisting in computing the ground state, then picking at random configurations 
and look at their distance to the ground state is not computationally practicable: configurations are far too numerous 
for this sampling to be efficient, let alone that one has to do it for numerous instances of the coupling constants in 
order to get relevant averaged quantities. We need to generate configurations whose energy density is somewhat under 
control. The technique we use consists in perturbing the Hamiltonian to increase the energy of the ground state, then 
compute the new ground state. The perturbation is chosen proportional to a small parameter e which we can tune 
to make the new ground state more or less distant from the original one: hereafter we will refer to this method as 
e-coupling. This is actually not a new idea 0,0,0,0,0, however to our knowledge so far this technique has 
been implemented only in numerical simulations. Here we show that an analytic solution is possible in the case of the 
simple random matching problem as well as the case of spin glasses on random graphs of fixed connectivity. 

In a recent paper Aldous and Percus [l5j used similar techniques to study the matching problem (both in the 
Euclidean and mean field versions) and the traveling salesman problem, conjecturing that it is possible to classify 
optimization problems into different universality classes according to the dependence of the ground state solution on 
small perturbations. We will compare our results to those presented in |l5| . and we will explain how the Replica 
Symmetry Breaking (RSB) transition can be detected within the cavity scheme. Let us point out that we will restrict 
our analysis to the level of the replica symmetric (RS) approximation although in principle there is no problem to 
extend the same analysis to higher level of replica symmetry breaking; we will discuss this possibility and the potential 
interest of this generalization in the context of other combinatorial optimization problems (e.g. SAT |2( and coloring 
i). 



2 



The rest of the paper has the following structure. In Sec. [H] we study the random simple matching problem: after 
a description of the model we show how the cavity approach works in general for this model and we introduce a 
generalization in order to deal with e-coupled systems. In Sec. IIIII we explain how the method works in the case of 
the diluted spin glass where a RSB transition is known to exist. In Sec. IIVI some conclusions and perspectives are 
presented. 

II. THE RANDOM SIMPLE MATCHING PROBLEM 
A. The Model 

Given an unoriented graph G = (V,E), where V are the vertices and E are the edges, a matching M is a set of 
edges having the property that no two edges in M have an end in common. We say that a vertex v £ V is matched 
if there is an edge incident to v in the matching. Otherwise the vertex is unmatched. A matching is called perfect if 
every vertex of G is matched. 

In the following we stick to the case where G is complete, i.e. there is an edge between any two vertices, and when 
we speak of a matching we mean a perfect matching. 

We introduce a function / defined on the set of edges E which associates a real number 1(e) to each edge e £ E. The 
quantity /(e) can be thought of as a distance, or - depending on the taste - as a weight, a cost, etc. We define the 
total length of a matching M as Lm = J2 e eM '( e )- The matching problem consists in finding the cheapest matching, 
i.e. which minimizes Lm- 

A classic example is the random Euclidean matching problem: the vertices are identified with N points drawn 
at random with the flat measure in the unit hypercube of a (/-dimension Euclidean space. The cost /(e) of an edge 
e = {v, w} is the usual Euclidean distance between the two vertices v and w. 

A mean field approximation of the Euclidean problem has been widely investigated [l6l fill ITsl Il9l I20L |2*H : the 
weights of the edges are i.i.d. random variables whose common probability distribution p is defined on the interval 
[0, +00 [. The function p(l) is assumed to behave like an integer power law for small /: 

p(l) ~ - (1) 
r! 

In the thermodynamic limit N — ► 00 limit, the mean distance of two nearest neighbors goes to like N~ s , where 
5 = l/(r + 1). Intuitively a minimum matching will only include edges of this order of magnitude, so that the only 
relevant feature of p is its behavior around 0, i.e. the r exponent. The one-edge and two-edge length distributions in 
this model match the ones of the Euclidean random matching problem in dimension d = r + 1 for short distances. 

Hereafter we will concentrate on the case r = 0. We will use the following conventions: there are N vertices; 
the distances lij = Iji between two vertices are i.i.d. random variables distributed following the flat distribution on 
[0,N] (this corresponds to a rescaling of a factor N with respect to Eq.QJ); as a shorthand notation, we indicate the 
set of these coupling constants as t. The length of the minimum matching is an extensive quantity, and the energy 
(Hamiltonian) of a matching can be defined equal to its length. A matching M can be univoqucly represented by a 
contact matrix (n,j) such that: 

• riij £ {0, 1} and fiy = n^. 

• nu = 0: no self linkage is allowed. 

• Vi, y^j—i n ij = 1- each site cannot be linked more than once. 

Obviously mj — if the edge {i,j} is not in M, and = 1 if the edge is in M. The entry is called the 

occupation number of the edge {i,j}. The Hamiltonian of the matching then reads: 

L t = n iJ l a> ( 2 ) 

l<i<j<JV 

where I plays the role of the quenched disorder, and the riy are the dynamic local variables. 

B. The cavity equations 

The cavity equations at finite temperature for the matching problem have been derived in l2l | (a comprehensive 
introduction to this subject is found in jisj)- Let us briefly reproduce the basic steps following |23j. The partition 
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function for the matching problem is 



z = ex p -w E 



Tlijlij 



(3) 



{Uij} 



l<i<j<N 



where the scaling factor N makes the free energy F = — \og[Z)/ j3 an extensive quantity. Following polymers theory 
|25|. one first introduces a more tractable representation of the partition function, which consists in mapping the 
matching problem onto a system of interacting spins. On each vertex i one puts a p-dimensional vector spin S, = 
(Si, . . . Sf) normalized by S^ = p. Let d\x be the integration measure on the corresponding sphere. If we define the 
coupling constants = exp (— /3iVZy), one can check that the partition function @ can be written as: 



Z = lim 



JV 



n^M(Si) 

.1=1 

Expanding the exponential into a power series and applying the following property: 



exp ( TijSi-Sj 

l<i<j<N 



lim / d /i (S i )5? 1 5?»..-S?« 

p^OJ 



(4) 



(5) 



one can easily recover (0 from |0J. Note also that the magnetization vector of spin i has components mf — # a ,i m i- 
The cavity method consists in adding a new spin So to a iV-site system {Si . . . Sat}. The new partition function is 
calculated assuming that the statistical correlations in the iV-site system can be neglected. More rigorously we make 
use of the clustering theorem, tacitly assuming that the system has just one pure state. We can thus encode the effect 
of the whole system onto each spin i as an effective field hi. The iV+l-site partition function can be written as: 



lim 



N 



(N N \ 

2J hiSi + r jS .Si J , 



(6) 



and can be easily computed thanks to giving: 



m 



JV 



(7) 



where we have used the fact that rrii, the magnetization of site i before the addition of site is 1/hi. At the end of 
the day we are not interested in the spin variables but in the solution of the matching problem. The thermal average 
of the occupation number of the edge 0— i is simply related to the magnetizations by: 



(not) = m T 0l m1. 



(8) 



Since we will be interested eventually in the ground state properties of this model we have to take the (3 — > oo limit 
of our equations. Following it is useful to set 

mi = e W< , for i G {1,...,N}. (9) 

The zero temperature limit of and |JSJ|, thanks to ©, reduces to the following zero-temperature cavity relation: 

(10) 

n 0i = , (11) 
where i* is the index attaining the minimum Nloi — fa. 



b Q = min Nl i - 4>i , 

i=l,...,N 



C. The e-coupling method 



The idea of the £-coupling method in the context of the matching problem is the following: given the set i of Zy 
distances (called 0-distances), one first finds the minimum matching (called 0-ground state), which is characterized 
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by some occupation numbers nij . Then one perturbs the lengths of the edges of the graph by adding a quantity e to 
the edges present in the 0-ground state. Formally the edge-lengths become the following e-distances: 



lij = kj +snij. (12) 

One solves the matching problem with these e-distances and obtains a solution we will name e-ground state, which 
is expected to be different from the 0-ground state. The larger e, the stronger the 0-ground state is penalized. The 
e-ground state is characterized by the occupation numbers nfj. Two quantities are of interest: first the difference of 
length (energy) between the e-ground state and the 0-ground state computed with the O-distances: 

A/., ^:,/,, n, ,)!,,. (13) 

i<j 

Second the distance d between the e-ground state and the 0-ground state: 

di = l- q e , (14) 

where the overlap qi is equal to the proportion of edges in common: 

2 

qe = ^J2 n ii n v- ( 15 ) 

i<j 

More precisely one would like to compute the average with respect to the coupling constants of these quantities, so 
let us define AL = ALg and q — qj. An analytic approach is possible in the framework of the cavity method. We 
have two spin systems, the 0-system and the e-system, standing on the same graph, and which are coupled. More 
precisely the Hamiltonian of the e-system is conditioned by the ground state of the 0-system. 

While the 0-system certainly obeys the cavity equations l|10l) . how to deal with the e-system is more problematic. 
A naive approach would lead us to: 

r = min 7i, (16) 

i— 1,...N 

where 

7i = hi -n if i ^ i*, (17) 
H=loi + f--Ti if i = i*. (18) 

Nevertheless these equations are wrong! They would be true if the interactions between the old spins in the e-system 
of N + 1 spins were the same as in the e-system of N spins, and this is not exactly the case: when adding the new spin 
to the 0-system, it gets matched to one of the old spins, whose previous match becomes unmatched. To circumvent 
this difficulty, one should distinguish between matched and unmatched spins in the 0-system. The variable Tj of the 
e-system will be called fa if the vertex i is matched in the 0-system, and Vi if it is not. The correct equations are thus: 

v Q = min (loi - /<), (19) 

2—1, ...N 

fo = . min 74, (20) 

I— 1,...N 

where 

Ji= hi-fi if M**, (21) 
ji=hi + e-Vi if i = i*. (22) 

The new spin in the e-system gets matched to the spin i**, the index which minimizes Eci. (|20|) . The contribution to 
the overlap q is <5j*,i**. The contribution to AL is hi** — hi*- 

When averaging over the disorder, in the thermodynamic limit the quantities (tpo, Vq, fo) and (<Pi,Vi,fi) are i.i.d. 
random variables (beware that tp, v and / on the same site are correlated). The above equations <|10[1 and l|19|) define 
a stochastic flow whose fixed point is the limit distribution of (ip,v,f). We use a population algorithm similar to 
the one discussed in 0,0 to solve the equations. In order to save computing time, we use a fluctuating connectivity 
approximation of the matching problem: we keep the only edges whose lengths are smaller than a given cutoff z, so 
that the connectivity of a vertex is a Poisson r.v. of mean z. We store a large population of M triplets (<pi, Vi, fi), i = 
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1, . . .AT, which we initialize randomly and update iteratively: at each step an integer k is extracted following the 
Poisson distribution of mean z, k random elements of the population are chosen and we compute (ipo, t>o, /o) following 
the scheme defined by Eqs. I(10() . (|19|1 and 1)20(1 limiting the search for the min to the k extracted triplets; the resulting 
(ipo, Vq, /o) overwrites an element of the population chosen at random. Once stochastic convergence of the population 
is achieved, we keep on iterating and compute en passant the contributions to AL and q. Their flat averages over 
many steps provide AL and d. The output of this algorithm is presented in Figs. Q an d J2J- Beside we have 
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FIG. 1: AL/N vs. d. The points obtained by the cavity population algorithm with z = 30, N = 200000, and 10 s iterations 
(CAVITY), coincide with the measures obtained by averaging over 4000 samples of total number of points N = 400 (TRUE). 

calculated the ground state of 4000 matching instances using BLOSSOM 4 software [2|| and we have tested the results 
against the cavity approach. In Fig. we display AL/N vs. d: the cavity approach is in perfect agreement with 
the direct calculation. In Fig. (J2J) we display the AL/N vs. d curve (main panel) and the best one-parameter fit 
of the form const x d 3 , while in the inset we display the d vs. e curve together with the best one-parameter fit 
of the form const x ^/e: in both cases the reduced x 2 is °f order 1. Note that a simple scaling argument shows 
that the two exponents are not independent: let us assume that for e small enough d ~ e@ , then from Eq. i|13|) 
AL/N ~ e l3+1 ~ d J {P+ 1 )/P . The scaling exponent a introduced by Aldous and Percus is easily recovered setting 
a = {j3 + l)//3. We find that the random matching features a scaling exponent f3 = 1/2 (a = 3), in agreement with 
the results found in A completely analytic study of the coupled system of Eqs. (|19fl and 120fl might be done, but 
we did not undertake it. 

An important remark is in order: in [24[ Aldous introduces and proves the following asymptotic essential uniqueness 
property: 

Let M be the generic element of a family of matchings depending both on N and on the realization of the Uj ; we call 
QM,min the overlap between M and the minimum matching, and dM,min = 1 ~ QM,min their distance. 
For each < S < 1 there exists e(5) > such that: if ViV, dm. min > S then 



(23) 



6 




In physical terms: in the thermodynamic limit, a matching (configuration) which differs from the ground state by 
a non zero proportion of edges has an energy density strictly greater than the one of the ground state. The other way 
round: a state with the same intensive energy density as the ground state can be obtained only by changing a non 
extensive number of edges in the ground state. This is the proof that there is no RSB, at least at zero temperature. 
The plot in Fig. can be seen as an illustration of this theorem. 

To conclude, note that our perturbation of the energy is O(N) so the information we get on the energy landscape 
is limited. In particular we do not explore the lowest lying excited configurations which have been numerically shown 
to have an energy AL ~ 1/V~N and d ~ 1/V~N |28j| . This is a limitation of our analytic approach: it is purely 
thermodynamic so that we cannot have a hint at finite size effects. 



III. THE GAUSSIAN SPIN GLASS ON THE BETHE LATTICE 



In close analogy with the previous section, we derive cavity equations for the £-coupling method applied to the 
Gaussian spin glass on the Bethe lattice of connectivity k + 1, in the presence of an external field. We keep to the 
level of the RS approximation. We will see that our results provide a self-consistency check of this hypothesis which 
enables to trace out when it is valid or not. 

Throughout this section we make thorough use of the notations introduced in our previous article We refer to 
sections II and III therein for details. 
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FIG. 3: A branch of a Cayley tree of connectivity k. 



A. The cavity equations for a single system at zero temperature 

First, following Mezard and Parisi [9j we work out the cavity method directly at zero temperature for a single 
system of N spins a interacting through coupling constants J - their set being denoted J - which are i.i.d. random 
variables with a Gaussian distribution of mean and variance 1, in the presence of an external field h ext . 

Let us consider the merging process of k branches rooted at spins i — 1, . . . k onto a new spin $ illustrated in 
Fig. ©. We look at how the energy of the ground state evolves under this process. Before the merging, on a given 
branch, we let the spin at the root undetermined. Thus the ground state energy of the branch rooted at spin i is 
conditioned by the value of the spin i and can be written: 

E(cr t ) = A l -h l a ll (24) 

where A{ is a constant, and hi is an effective field (beware that it is not the local field). Note that hi contains the 
effect of the external field h ext . As a consequence, the energy of the system of the k branches before the merging is 
conditioned on the values of the spins 1, . . . k, and reads 

E(a u . . . <T fc ) = Ax + . . . + A k - hiai - ... - h k a k . (25) 

The system after the merging has an energy 

E'(ai, . . . a k ,cr< s >) = Ai + . . .+A k -(hi + Jq >t ia^)a 1 -.. .-(h k + J^ ik aq,)a k ~h ext cr^, (26) 

which, in order to be that of the ground state, is to be minimized with respect to a\, . . . a k at fixed cr$. This is realized 
by independently choosing the sign of each ai such that {hi + J^,^a^)ai = \h% + Jq,,i<j§\. As we can write 

\hi + J*,iO$| = Lo(hi, J<s> ti ) + X(hi, J*^) "*) ( 27 ) 

where 

„ M = fcM+fcil, (28) 



we get the appealing following form: 



k I k \ 

cr$ . (29) 



By comparison with Eq. (|24|l this leads to the recursion relation: 



B. The cavity equations for two coupled systems 
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FIG. 4: q vs. e for k — 2, in the case h ext = 0. 

Now we imagine there are two spin systems sitting on the N vertices A of the same Bethe lattice: 

• the a spins obey the same Hamiltonian as in the previous section Tin Al 

H [a] — — Ja,b&a&b ~ h ext a a- 

(A,B) A 

Let us call a* — the ground state of this Hamiltonian. 

• the r spins obey the following perturbed Hamiltonian, conditioned on the <j* a : 

H"[t|<7*] = - Ja,btatb ~ K xt ta -e^2<T* A T A . 

A,B A A 

The choice e > corresponds to an attractive interaction for the r spin variables to the configuration a* 
choice e < to a repulsion. The ground state of the Hamiltonian (|32|) is named r* = (t a ). 
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FIG. 5: q vs. e for k = 2, in the case h ext = 0.1, obtained with Af = 2000 and 100000 x Af iterations. 
We are interested in the overlap 



q = qj = — V't. 



1 1 



(33) 



where with the symbol Oj we mean the average of a generic j7-dependent observable O j over the different sample 
realizations. The correlations between the a spins and the r spins can be accounted for by stating that the value of 
the effective field acting onto spin t$ depends on the value of spin c$ , where $ is the root of an isolated branch. So 
we have to store three quantities related to the root $ of a branch: 

- h®, the effective field acting onto spin cr$; 

- /ij, the effective field acting onto spin r$ under the condition that spin cr$ has the value +1; 

- the effective field acting onto spin r$ under the condition that spin cr$ has the value -1. 

1. Iteration 

When performing the merging of k branches onto a new site, the equation for the system of the a spins is exactly 
the same as in the previous section UlI Al 



h$, = ^2 X(hi, J$,i) + h e 



(34) 



i=l 
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FIG. 6: q vs. h ext for k = 2, e = 0, obtained with N = 2000 and 100000 x N iterations. 



As far as the r spins are concerned, one must compute /ij as well as h^. So we assume in turn that <r$ = +1 and 
ct$ = — 1 . This determines the values of the a* : 

a* = sign(/i. t + J$,i(T$), (35) 

and the effective field acting onto spin t$ reads 

k 

h °i = J2 x ( h ? > J i) + £(J * + h ^t- (36) 

i=l 

2. Measures 

To measure the overlap, one uses the merging procedure of (k + 1) branches onto a new vertex the effective field 
acting onto spin aq, is: 

fc+i 

Hq, = '^2\(h il Jq, ti ) + h ext . (37) 

i=l 

By contrast with the iteration, spin a^j has a determined value, the one which minimizes the energy: erj =sign(.ff^). 
This also determines the values of the a* according to Eq. Il-ioll . Eventually one can compute the effective field acting 
onto spin t^: 

k+l 

H* = J2 KK l , Ji) + + hext, (38) 
i=i 

whose sign gives t£. The contribution to the overlap is ctJt^. 
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C. The results 

Unless otherwise stated the following results are for k = 2. Once again they have been obtained by a population 
algorithm. Here the population is made of M triplets (h, h + ,h~). 

In the case e > we expect q to be equal to 1, regardless of the value of h ext : the ground state of the perturbed 
Hamiltonian (|32[l should be the same as the Hamiltonian l|31|) . We will see that so it is. Thus the interesting regime 
is e < 0. 

What we obtain in the case h ex t = is plotted in Fig.Q. The fact that q = — 1 for every e < is due to 
the symmetry of the original Hamiltonian (|31|) under reversal of all the spins: the ground state of the perturbed 
Hamiltonian is t a = a* A . More interesting is the fact that at e = 0, q = 0: assuming RS one would expect q — 1. This 
is a sign that the RS ansatz is not self consistent, and is to be dismissed. 

Things may be more convincing if one turns the external field on. The case h ex t — 0.1 is plotted on Fig. 101. As the 
symmetry of the original Hamiltonian is lifted, the plot of q as a function of e < is no longer a constant equal to -1. 
It does however tend to -1 when e — > — oo because in this limit the attracting term between the a system and the r 
system dominates the perturbed Hamiltonian l|32l) . When e — ► 0~ q goes continuously to a value which is no longer 
0, but is still not 1. Again the RS ansatz is in trouble. 

We do expect however the RS ansatz to be valid for a sufficiently high external magnetic field. Fixing e — 0, we let 
h ext increase. See Fig.©: it appears that q is an increasing function of h ex t and saturates to 1 at h c ext ~ 0.48. It is 
the sign that the RS ansatz becomes self-consistent above h% xt . 

For fc = 5we find h c ext ~ 1.86 which is different from the value h c ext ~ 2.1 found in 30] by the analysis of numerical 
simulations. We believe our result is exact, and the discrepancy can be explained by the fact that their result relies 
on finite size scaling arguments with a relatively poor precision. 

What we got here is actually the point at T = of the de Almeida-Thouless (AT) line [3l| for the Gaussian spin 
glass on the Bcthc lattice. We had the idea we could generalize our approach to a non zero temperature to determine 
the whole line. 



D. Searching for the AT line 

It is not difficult to generalize the argument of section IIII Bl fin the only case e = 0) to a non zero temperature. We 
consider two non interacting systems, the a spins and the r spins, standing on the same Bcthc lattice and obeying 
the same Hamiltonian l|31l) . Solving this double system at the level of RS is easy: we use the population algorithm 
of section IIII Al adapted to follow simultaneously two populations. A crucial point is that whenever one randomly 
extracts sites or coupling constants, they are the same for the a population and the r population. This procedure 
enables to measure the average overlap between the two systems 

q = qj = (<t a )(t a ). (39) 

The criterion for the RS ansatz to be self-consistent is 

q = ^5, (40) 

where m A is the local magnetization measured either in the cr-system or the r-system (they are obviously equal). 

Given a value of T, we run the algorithm for increasing values of h ex t, so as to determine its value h c ext {T) beyond 
which condition (gUJ holds. The plot h c ext (T) is the AT line. See Fig.Q). To our knowledge it is the first time the AT 
line has been obtained for a spin glass on the Bethe lattice. Two predictions made in [32j can be checked. First, the 

critical temperature T c , such that h^ xt (T c ) = 0, is the solution of the equation k tanh 2 ( J/T c ) = 1: in the case k = 2 
this gives T c = 0.748. Second close to T c , h c ext (T) should behave like (T c — T) 3 ' 2 , which a numerical fit of our data 
confirms. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this paper we presented the derivation and implementation of the e-coupling method in the framework of the 
cavity approach. This technique allows us to explore the way the energy of a configuration varies as a function of 
its overlap with the ground state and more generally to address the problem of the organization of the lowest energy 
configurations. The lesson we get applying this method to the case of the simple random matching problem is very 
clear: the space of the lowest energy configurations is organized such that their energy difference with respect to their 
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FIG. 7: The de Almeida-Thouless line (h vs. T) on the Bethe lattice in the case k — 2, separating Spin-Glass (SG) phase from 
the Paramagnetic (PM) phase. Data are obtained with N = 10000 and 10000 x N iterations. We estimate error bars to be 
hardly visible on this scale. Continuous line is just a guide to the eye. 



distance from the ground state scales as AE/N oc d 3 . A situation like this, or in general whenever AE/N oc d a 
with a > 0, is related to the property of replica symmetry of the system, which implies Aldous asymptotic essential 
uniqueness property |24j . 

A similar study presented in |T^ suggest that this property is also shared by the Minimum Spanning Tree problem, 
the Minimum Matching Problem in Euclidean dimension d — 1, and the Traveling Salesman Problem also in Euclidean 
dimension d — 1 (all with a — 2). Minimum Matching Problem and Traveling Salesman Problem in d = 2,3 are 
instead characterized by a — 3 as the mean field matching problem we have studied. 

A simple case with no asymptotic essential uniqueness property is the spin glass on a fixed connectivity random 
graph studied in Sec. IIHI Indeed our computation based on the RS assumption yields that lim e _^ - d ^ which is 
a physical nonsense since the model has a unique ground state (the couplings are gaussian). This inconsistency tells 
us not only that - as we already know after [3, 13 - the cavity approximation must be improved in order to take into 
account the presence of many states, but also gives us a practical tool to probe the phase space for the onset of full 
RSB: the search for the AT line in the case of the spin glass on a fixed connectivity random graph, is a simple and 
instructive example. 

A very interesting issue is the generalization of the e-coupling method to the case where the asymptotic essential 
uniqueness does not hold. In the last year a compact and efficient formalism has been developed to apply the cavity 
method to SAT and Coloring Problems 0,0 , at the level of 1RSB, evidentiating a clustering transition. This clustering 
transition consists in the sudden appearance of an exponential number of metastable states, which - intuitively - cause 
local search algorithm to get stuck. We believe that it is possible to generalize the £-coupling method presented in 
this paper to the 1RSB level, although at a higher computational cost. This could give interesting results on the 
inner mechanism of the clustering of states. Open problems like how large a single cluster could be and what the 
inter-cluster mean distance is, could be addressed within this formalism. 
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